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Abstract 



Aims. We seek to understand the propagation mechanisms of positrons in the interstellar medium (ISM). This understanding is a key 
to determine whether the spatial distribution of the annihilation emission observed in our Galaxy reflects the spatial distribution of 
positron sources and, therefore, makes it possible to place constraints on the origin of positrons. 

Methods. We review the different processes that are likely to affect the transport of positrons in the ISM. These processes fall into 
three broad categories : scattering off magnetohydrodynamic waves, collisions with particles of the interstellar gas and advection with 
large-scale fluid motions. We assess the efficiency of each process and describe its impact on the propagation of positrons. We also 
develop a model of positron propagation, based on Monte-Carlo simulations, which enable us to estimate the distances traveled by 
positrons in the different phases of the ISM. 

Results. We find that low-energy (< 10 MeV) positrons generally have negligible interactions with magnetohydrodynamic waves, 
insofar as these waves are heavily damped. Positron propagation is mainly controlled by collisions with gas particles. Under these 
circumstances, positrons can travel huge distances (up to ~ 30 kpc/;7jj^^,-3 for 1 MeV positrons) along magnetic field lines before 
annihilating. 



1. Introduction 

Positron annihilation in the Galactic center (GC) region is 
now a firmly estabUshed source of radiation, which has been 
observed since the early seventies in several balloon and satel- 
lite experiments ( see [win Ballmoos et al. 2003, Jean et al. 2004, 
iDiehl et al.l2006l for reviews). Despite significant progress in ob- 
servational capabilities, the origin of Galactic positrons remains 
an open question. 

Recent observations of the 511 keV line intensity using 
the SPI spectrometer onboard the International Gamma-Ray 
Laboratory (INTEGRAL) observatory have revealed a dif- 
fuse em ission, distributed in the bulge and the disk of ou r 
' Galaxy dKnodlseder et al.l 120051: IWeidenspointner et al. I I2OO6I) . 
Observations further indicate that the bulge-to-disk luminosity 
ratio (hereafter B/D) of the 511 keV line is rather large 3-9) 
compared to the distribution of any candidate source. Under the 
hypothesis that positrons annihilate close to their sources, the 
spatial distribution of the annihilation emission should reflect 
the spatial distribution of the positron sources. In this view, the 
large B/D ratio could be explained by sources belonging to the 
old stellar population dKnodlseder et al.ll2005h and a disk emis- 
sion, partly or totally, attributed to the radioactive decay of -^Al 
and '^'^Ti produced in massive stars. The most recent SPI ana- 
lysis covering more than 4 years of data shows hints of a lon- 
gitudinal asymmetry in the spatial distribution of the 511 keV 
line emission produced in th e inner part of the Galactic disk 
dWeidenspointner et al.ll2008l) . A similar asymmetry is observed 
in the distribution of low-mass X-ray binaries emitting at high 
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energies, suggesting that these objects might be the dominant 
sources of positrons. This conclusion is contingent upon the hy- 
pothesis that positrons annihilate close to their sources, a hypo- 
thesis that should be called into question. 

Several authors argued that positrons annihilate in the vici- 
nity of their production sites, based on the assumption that par- 
ticles propagate according to the so-called Bohm diffusion, i.e., 
with a mean free-p ath equal to their Larmor radius dBoehm et al.l 
llOOllmngllOOel) . This implicitly supposes that magnetic fluc- 
tuations in the ISM are strongly tangled on all scales, as is li- 
kely the case near strong shocks in supernova remnants. In rea- 
lity, given the complexity and the variety of thermodynamical 
phases in the ISM, Bohm diffusion pro b ably overestimate s par- 
ticle confinement dParizot et al.ll2004al) . iJean et all d2006l) esti- 
mated the distances traveled by positrons using a model combi- 
ning a quasi-linear diffusion theory of wave-particle interactions 
at high energy and an approximate propagation model including 
collisions with the ambient ISM matter at low energy. A prelimi- 
nary model of the gas content in the Galactic bulge (a more com- 
plete descri ption of the gas spat ial distribution in this region is 
presented in lFerriere et al.ll2b07h was then added to the propaga- 
tion model. The authors concluded that MeV positrons injected 
by radioactive processes into the bulge cannot escape from it and 
that a single source releasing positrons in the bulge might have 
difficulties accou nting for the obse rved spatial extent of the anni- 
hilation emission. |Prantzosld2006l) proposed a solution to explain 
the large B/D ratio measured with SPI, in which positrons pro- 
duced by type la supernovae in the old disk are tr a nspor ted along 
magnetic field lines into the bulge. ICheng et all d2006l) sugges- 
ted that positrons in the bulge could originate from the decay 
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of produced in high-energy pp colUsions. In this work, the 2.1. Positron resonance with MHD waves 



energetic protons would be shock-accelerated when a star is ti- 
dally disrupted in the accretion disk of Sagittarus A* (Sgr A*). 
Since high-energy (> 30 MeV) positrons take as long as ~ 10^ 
years to cool down, they should be able to propagate far away 
from Sgr A* and to fill the Galactic b ulge. In order to explain 
the annihilation emission in the bulge, iTotanil ( |2006|) proposed 
a scenario in which MeV positrons were produced in the accre- 
tion disk of Sgr A* ~ 10^ years ago. These positrons would have 
filled the bulge while being transported by large-scale outflows 
(~ 100 km s"'). All these studies invoked positron transport to 
explain the observed annihilation emission, but they did not in- 
clude any detailed examination of the physical processes that 
could contribute to it. 



More recentlv. lHigdon et al.l (|2009|) discussed the role of ma- 
gnetohydrodynamical (MHD) fluctuations in positron transport 
in the different ISM phases. They further estimated the positron 
mean free-path using a phenomenological model of MeV elec- 
tron transport in interplanetary turbulence, especially in the ioni- 
zed phases where turbulence is undamped down to small scales 
(see appendix). However, the model used to derive the positron 
mean free-path is strictly valid for interplanetary plasmas, and 
the nature of ISM turbul ence is likely different. The model de- 
veloped by iHigdon et al.l ( i2009i) is up to now the most valuable 
attempt to account for the propagation of low-energy positrons 
in the ISM. Here, rather than considering a particular turbulent 
model, we review the transport processes that govern the posi- 
tron mean free-path. 

We discuss three different regimes of positron transport in 
the Galaxy, namely, the regime dominated by scattering off 
MHD waves (section|2]), the regime dominated by collisions with 
particles of the interstellar gas (section [3]l and a regime of ad- 
vection with large-scale fluid motions (section |4|i. Based on the 
quantitative results obtained in each section, we infer the rele- 
vance and the importance of the respective regimes. In section|5j 
we summarize our study and discuss its possible implications for 
positron transport at Galactic scales. 



2. Scattering off magnetohydrodynamic waves 



In section 12. H we present the condition for wave-particle 
resonance, restricting our analysis to magnetohydrodynamic 
(MHD) waves. We further discuss the effective value of the 
Alfven speed in the neutral phases of the ISM. In section IZ21 
we present the properties of the MHD wave cascades initiated 
at large spatial scales, e.g., by the differential rotation of the 
Galaxy or by the explosion of supernovae. In section 12.2.11 we 
present the different energy transfer times associated with the 
different type of cascades that likely developped in the ISM. In 
section [2.2. 21 we discuss the collisional and collisionless dam- 
ping mechanisms of MHD waves as well as the smallest scales 
of the MHD cascades in the different ISM phases. Only a brief 
summary is given there ; the technical derivation is postponed to 
AppendixlAl In section [231 we discuss the global effect on posi- 
tron transport of MHD waves injected at large scales. In section 
12.41 we consider the local effect of magnetic fluctuations gene- 
rated by the streaming instability, which can survive down to the 
small scales of low-energy positrons. 

Throughout section[2] we assume, for simphcity, that the in- 
terstellar gas contains only hydrogen. 



In a medium with non-vanishing magnetic field, wave- 
particle interactions proceed through the Landau-synchrotron re- 
sonance condition expressed as (Melrose 1986.) : 



(1) 



where co is the wave frequency, the component of the wave 
vector along magnetic field lines, the particle velocity along 
field lines, £ an integer, and Qse = i^ce/y the synchrotron 
frequency of the particle (a positron e^ in the case at hand, 
hence the subscript e), with y the particle Lorentz factor, Oce = 
{eB/mgc) the particle cyclotron frequency, and B the magnetic 
field strength. 

In this work, we consider only MHD waves, more specifi- 
cally, shear Alfven waves and fast magnetosonic waves as these 
waves are a m ajor part of the magnetic flu ctuations that per- 
vade the ISM dLifliwick & GoldreichlllOOl . We do not consi- 
der slow waves separately, because the dynamics of the slow 
wave cascade were shown to be entirely controlled by the Alfven 
wave cascade dLithwick & G oldreich 2001) and the slow wave 
spectrum is basically the s ame a s the Alfven wave spectrum. 
Following lYan & LazarianI (|2004|) . we assume that the transport 
of charged particles is governed by either Alfven waves or fast 
magnetosonic waves. Note that higher-frequency waves are po- 
tentially important as well, as they can easily fulfill the above 
resonance condition with positrons produced by radioactive de- 
cay (energy < 1 MeV). However, whistler waves, the most in- 
teresting waves in this frequency domain, are right-handed po- 
larized and, therefore, cannot be in resonance with positrons. 
Higdon et al. (2009) have discussed the possible role of self- 
generated whistler waves. We will come back on this possibility 
in section 12.41 The possible effects of large-scale compressible 
motions will be examined in section|4] 

For the waves of interest here, it can be shown that co «: 
QseE so we are entitled to neglect o) in Eq. ([T]i, except in the 
particular case { = 0. The case £ = corresponds to the so- 
called Cherenkov resonance, in which the particle interacts with 
a zero-frequency wave in a frame moving at velocity v^\. In that 
case, Eq. ([T) reduces to w = The Cherenkov resonance can 
be important for magnetosonic waves, which have a perturbed 
magnetic field component parallel to the background magnetic 
field. 

Here, we consider only the dominant harmonics i = +1. 
With cl) «: i^se, Eq. ([T]i can then be recast into the form 



k\\ r'^ cos a ^ ±\ 



(2) 



where a is the particle pitch angle (angle between the particle 
velocity and the local magnetic field), r'^ - {vjO.^^) - (pc/eB) is 
the particle gyroradius divided by sin a (simply referred to as the 
gyroradius in the following), and p = ynif^u is the particle mo- 
mentum. Eq. (|2]i expresses the resonance condition in terms of 
the particle gyroradius. It can be rewritten in terms of the particle 
momentum as 

eB 

P = ±-r (3) 

ck\\ cos a 

or in terms of the particle kinetic energy, Ek, as 



Ek = (me c ) 



eB 



I nieC- k\\ cos a 



+ 1-1 



(4) 



' According to Eq. this condition is automatically satisfied for 
positrons with y « nip/in^, corresponding to energies « 0.94 GeV. 
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Roughly speaking, MHD waves can exist only at frequencies 
lower than the proton cyclotron frequency, Qcp - {eB/nipc) : 



CO < Q.r 



(5) 



Either they are damped by collisional effects (mainly viscous 
friction and ion-neutral collisions) at low frequencies (see sec- 
tion lA.lb or, if they manage to survive collisional effects, then 
at frequencies approaching Ocp, they are heavily damped by 
Landau damping due to thermal protons (see section [ASI i. 

For Alfven waves {co - VAk\\, with Va the Alfven speed), 
Eq. (|5]l is equivalent to 



^11 < 



■cp 



(6) 



which, in view of the resonance conditions, Eqs. (O and (HI, 
implies a threshold on the positron momentum (obtained for 
I cos a| = 1) : 

p > mp Va , (7) 
and a threshold on the positron kinetic energy : 



Ek > (OTeC^) 



Me C 



+ 1-1 



(8) 



For fast magnetosonic waves (w - V^iQ) k, with the angle 
between the wave vector and the magnetic field, Vp the phase 
speed of the fast mode, comprised between Vf(0) - max(yA, Cs) 

and V'F(;r/2) = V^s = ^V^^Tc^, and Cs the adiabatic sound 

speed), Eq. (|5]l leads to ^ < Qcp/Vf(0) and, hence, p > 
mp yF(0)/ cos 0. Since the functions yF(0) and yF(0)/cos0 
reach their minimum values at = 0, the absolute requirement 
on the wavenumber reads 



k< 



■cp 



Vf(0) ' 



(9) 



and the thresholds on the positron momentum and kinetic energy 
(again obtained for | cos ffl = 1) are given by 



p > mp yF(0) 



and 



Ek > (WeCO 



f mp yF(0) 



1 - 1 



(10) 



(11) 



respectively. 



The question now is what expression should be used for the 
Alfven speed. In a fully ionized medium, the Alfven speed is 
simply Va - B/ -sJ^Kpi, with pi the ion mass density. However, 
in a partially ionized medium, the relevant Alfven speed de- 
pends on the degree of coupling between ions and neutrals. If 
the ion-neutral and neutral-ion collision frequencies, Vin and Vni, 
are much greater than the wave frequency, ai, then ions and neu- 
trals are very well coupled together through ion-neutral colli- 
sions, and as a result, an Alfven wave will set the entire fluid 
(ions + neutrals) into motion. In that case, one should use the to- 
tal Alfven speed, yA.tot = -^J^^Ptou with = pi + Pn the total 
(ion + neutral) mass density. In contrast, if the ion-neutral and 
neutral-ion collision frequencies are much smaller than the wave 



frequency, then ions and neutrals are no longer coupled together, 
and an Alfven wave will solely set the ions into motion. In that 
case, one should use the ionic Alfven speed, yA.i = B/ ^j^^p■^, as 
in a fully ionized medium. 

In an atomic medium with temperature T < 100 K, the 
ion-neutral and neutral-ion collision frequencies are given by 
Vin - (1-6 X 10"^ cm-* s"')«n and ^ (1.6 x lO"** cm-' s"')ni, 
where and rij are the neu tral and ion number densities, res- 
pectively jOsterbrockll 1 96 Ih . At high temperature, the collision 
frequencies increase as Vf (lBraginskiilll"96l IShuU & Draing 
1987.) ; assumi ng an effective cross section for H-H^ collisions 
~ 10"^ cm^ (IWentzelllT97l . we find that, for T > 140 K, 
Vin ~ (1-4 X 1 0''' cm^ s'') Vr/100K«n and Vni ~ (1.4 x 
10'^ cm"* s"') V 7^/100 Kni. In a molecular medium, the colli- 
sion frequencies are Vjn - (2.1 x 10~^ c m^* s"')/in and Vni - 
(2.1 X 10 cm^ s ')ni (lOsterbrockll 1961b . For comparison, the 
frequency of a resonant Alfven wave is given hy u - Va^w 
with, according to Eq. (|2]i, k\\ > l/r'^ - (eB/pc), i.e., o) > 

(6.6 X 10"^ s"') (fi^c/ V"cin-^) iP^^MeV whcrc n is the relevant 
number density (e.g., n - rti + n„ if Vm, Vni » w and n - tiiif 
^in,Vni <y). From this, it follows that for typical interstellar 
conditions (see Table [U, Vin, Vni <«; w up to positron energies of 
at least 1 GeV. Hence, in the present context, the relevant Alfven 
speed is yA,i - B/ ^j'\np■^, not only in the ionized phases, but also 
in the so-called neutral (i.e., atomic and molecular) phases. 

Alfven waves 

With the above statements in mind, the requirement on the 
parallel wavenumber of Alfven waves, Eq. (l6|, can be rewritten 

as 



h < 



= (4.4 X 10"** cm"') 



(12) 



The corresponding condition on the positron kinetic energy, 
Eq. ([8]l, becomes 



Ek > (511 keV) 



|l.8x 10-4 + 1-1 



(13) 



The above expressions were obtained on the assumption that the 
only ion present in the ISM is H^. To account for the presence 
of other ions, it suflices, in good approximation, to replace «i 
by + 4 hhc with «He -0.1 hh in t he hot phase and by «h+ + 
12 «c,gas with ncgas - 1.4 X lO ^* nil dCardelli etalJI 19961) in the 
atomic phases. 

The maximum parallel wavenumber, fcymax, and the mi- 
nimum kinetic energy, isk.min, given by the right-hand sides 
of Eqs. ( fT2] i and (fTsT l. respectively, are listed in Table [T] for 
the different phases of the ISM. Also listed in Table [T] are the 
estimated temperature, T, magnetic field strength, B, hydrogen 
density, riu, and ionization fraction, /ion - «i/(«i + Hn) of 
the dififerent phases. The values of T, hh and /ion are taken 
from the review paper of iFerrier i (I2001h . For B, we adopt 
the value of 5 infer r ed from rotation measu re studies (e.g. 
iRand & Kulkarnil Il989t lOhno & Shibatal Il993h for the warm 
phases, t he value of 6 nG iriferred from Zeeman splitting measu- 
rements (iHeiles & Trolandl2005b for the cold phase, th e relation 
B PC ^Jnu normalized to B = 85 i-iG at tin -10'^ cnT^ dCrutcheJ 
1 19991) for molecular clouds, and the two extreme values of 
2 fiG and 20 pG for the hot phase. The lower value pertains 
to the standard scenario in which the hot gas is generated by 
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5 


0.2-0.5 


0.007 - 0.05 
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50 - 100 


6 
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0.008 - 0.05 
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MM 


10-20 
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10^ - 10*^ 
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Tab. 1. Physical parameters of the different ISM phases : molecular medium (MM), cold neutral medium (CNM), warm neutral 
medium (WNM), warm ionized medium (WIM) and hot ionized medium (HIM). T is the temperature, B the magnetic field strength, 
«H the hydrogen density, /ion - «i/(ni + nn) the ionization fraction, rii the ion density, A;||niax the maximum parallel wavenumber of 
Alfven waves (right-hand side of Eq. (fT2]i). and £'k,min the minimum kinetic energy required for positrons to interact resonantly with 
Alfven waves (right-hand side of Eq. (fT3]l). Here, we assume a pure-hydrogen gas, for which ni = /ion hh and «n = (1 - /on) "h- 



stellar winds and supernova explosions, which sweep up the 
ambient magnetic field lines and evacuate them from the hot 
cavities. The higher value pertains to an alternative scenario in 
which large-scale highly turbulent MHD fluctuations produce 
magnetic fields above equipartition with th e local thermal 
pressure (iBykovl 1200 ll : iParizot et al.l l2004bl : and references 
therein). 

From Table [1] it emerges that the maximum parallel wave- 
number of Alfven waves is typically a few 10"'' cm"', close 
to the largest wavenumber, k]^ ^ 10"' cm ', of the elec- 
tron density power spectrum inferred from interstellar scintilla- 
tion ([Armstrong et al.l ll995^. Furthermore, the minimum kinetic 
energy required for positrons to interact resonantly with Alfven 
waves varies from a few keV (in the warm ionized medium) to 
a few hundreds of keV (in regions with large Alfven speeds, na- 
mely, in molecular clouds and possibly in the warm neutral and 
hot ionized media). For comparison, positrons produced by ra- 
dioactive decay are injected into the ISM with typical kinetic 
energies ~ 1 MeV. This means that positrons from radioactive 
decay can interact resonantly with Alfven waves only over a 
restricted energy range. This range is particularly narrow in re- 
gions with large Alfven speeds, such as molecular clouds (where 
Va,! S 185 kms"') ; it can even vanish in the hot phase if the ma- 
gnetic field is as strong as 20 yuG (implying VA,i > 436 kms '). 
In contrast, the resonant range extends over at least two orders 
of magnitude in the warm ionized medium, where the ion den- 
sity is highest (and VA,i - (16 - 31) kms"'). The reason why 
resonant interactions with Alfven waves are no longer possible 
below /ik.min IS becausc the Larmor radius has become smaller 
than the smallest possible scale of existing Alfven waves. 

Fast magnetosonic waves 

The numerical expressions and values of the wavenumber 
and kinetic energy thresholds can be obtained in the same man- 
ner as for Alfven waves, with the two following differences : 
First, the wavenumber threshold (Eq. ^) applies to the total 
wavenumber, as opposed to the parallel wavenumber Second, 
the speed entering the expressions of the thresholds is the phase 
speed of the fast mode for parallel propagation, V'f(O), as oppo- 
sed to the Alfven speed, Va- In practice, however, the second 
difference is only formal, except in the hot low-B phase. Indeed, 
Vf{0) = max(VA, Cs), and max(yA, C,) = Va in all the ISM 
phases, except in the hot low-B phase, where max(yA, Cs) - 
Cs - 166 kms '. In consequence, the values of the maximum 
wavenumber, kn^^x, and the minimum kinetic energy, fik.min, are 
those listed in Table [T] except in the hot low-B phase, for 
which Eqs. (|9]l and (fTTT i lead to fc,nax - 1-2 x lO"** cm"' and 



Ek,min - 220 keV, respectively. The latter value is roughly an 
order of magnitude larger than for Alfven waves, which signifi- 
cantly narrows down the energy range over which positrons from 
radioactive decay can interact resonantly with fast magnetosonic 
waves. 

The Cherenkov resonance (for f - 0) occurs when co = 
k^ uy, i.e., at a wave propagation angle ©t such that Vp{@j) = 
COS0T, independent of the wavenumber. Since Vp{®t:) ne- 
ver departs from the fast magnetosonic speed, Vms, by more 
than a factor V2, this expression is approximately equivalent to 
COS0T - Vms/y\\- The Cherenkov resonance requires that the 
wave at propagation angle ©t exist and not be damped by a col- 
lisional or collisionless process (see next subsections). 

Let us re-emphasize that the above results should be consi- 
dered only as rough estimates. As we will now see, both Alfven 
and fast magnetosonic waves are subject to various damping pro- 
cesses in the ISM. 

2.2. MHD wave cascades 

2.2.1. Energy transfer timescales 

We consider both Alfven waves and fast magnetosonic 
waves as parts of turbulent cascades. The main sources of tur- 
bulence able to counterbalance the dissipation mechanisms ex- 
pected in the ISM are the magnetorotational instability driven 
by the differ ential rotation of the Galax y and the explosion of 
supernovae 

dM ac Low & Klessenll2004l) . Both mechanisms re- 
lease energy at large scales. In the rest of the paper, we adopt 
Linj = 100 pc for the injection scale. This does not preclude the 
possibility of injecting magnetic fluctuations at smaller scales ; 
an issue discussed in section [241 

Alfven wave cascade 

The most recent developments in MHD turbulence theory 
explain the energy cascade towards smaller scales by the dis- 
tortion of oppositely trave lling Alfven wave packets (e.g., 
iLithwick & Goldreichll200lb . The kinematics of the interactions 
produce a highly anisotropic cascade, which redistributes most 
of the energy in the perpendicular scales^ We will return to this 
important question in section [273] 

The transfer time of the Alfven wave cascade, ta, corres- 
ponds to the wave packet crossing time along the mean magnetic 



^ Throughout this paper, the perpendicular and parallel directions are 
taken with respect to the mean magnetic field direction. 
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field: 



For reference, when Linj = 100 pc, one has 



ta = (1.4 X 10'^^ s) 



A: Linj cos© 



(14) 



where n is the relevant number density (in a fully ionized me- 
dium, n - Hi, while in a partially ionized medium, n - ni + tin if 

Vin, Vni » OJ and « = «i if Vin, Vni « Oj). 



for the Alfven cascade and 



Ttf = 1 



(18) 



for the fast magnetosonic cascade. If the collisional damping 
rate, Fcoii, is high enough that Fcou ta = 1 at a scale larger than 
/Ip, then the Alfven cascade is cut off by collisional damping at 
that scale. On the other hand, if Fcoii ta < 1 down to A^, then 
the Alfven cascade proceeds down to the collisionless range and 
is eventually cut off by collisionless damping at a scale smaller 
than /Ip. Similarly for the fast magnetosonic cascade. 



Fast magnetosonic wave cascade 

The transfer time of the fast magnetosonic wave cascade, rp, 
can be written as 



Tf = 



OJ 



1/2 
inj 



6v^. 

inj 



-1/2 



(15) 



where 6vk is the turbulent velocity at scale A:"', dviai is the tur- 
bulent velocity at the injection scal e, Ljnj, and, as before , Vp 
is the phase speed of the fast mode dYan & Lazarianll2004l) . In 
writing the second identity, we assumed a Kraichnan spectrum 
(dE/dk oc k^^/~ or Svi^ oc A:"'/"*). For our numerical estimates, we 
adopt again Linj - 100 pc and we set both Svi^j and Vp to the fast 

magnetosonic speed, V^s = -^V^ + Cl (as mentioned earUer, Vp 

does not depart from V^s by more than a factor V2), whereupon 
we find for L; 



■inj 



100 pc 
Tp = (1.4 X 10'^ s) 



(^^inj) ^ 



with the density n defined as for the Alfven wave cascade. 

For both Alfven and fast magnetosonic waves, the dominant 
damping process depends on the wavelength (or inverse wave- 
number) compared to the proton collisional mean free-path. 



j2 

VpTp ^ (3.5x 10'^ cm) — — 



Hp 



A 



(16) 



where Vp - -y/T^BTpTwip is the proton r.m.s. velocity, Tp the 
proton collision time, «e the electron density, and A the Coulomb 
logarithm, given by 



and 



A 



25.3 - 0.5 Innecm -i + In T'e.ev 



for Te < 50 eV 



for Le > 50 eV 



(lBraginskiilll965h . Waves with k ^ > Ap are basically collisional 
and, therefore, affected by collisional damping (see section lATTT i. 
whereas waves with A: ' < Ap are basically collisionless and af- 
fected by collisionless damping (see section |A!2] i. 

As we will see below, in all the cases considered here. 



-inj 



> Ap, which means that both turbulent cascades start in the 



collisional range. The wave energy is then transferred to smaller 
scales up to the point where the wave damping rate, F, becomes 
equal to the transfer rate. In other words, the turbulent cascades 
are cut off at a wavenumber k^ut such that 



Fta = 1 



(17) 



2.2.2. Damping and cutoff 

In this section, we summarize the investigation presented 
in Appendix |A] on the dominant damping processes of the 
Alfven and fast magnetosonic wave cascades in the different 
ISM phases. 

In the mostly neutral, atomic and molecular phases of the 
ISM, the Alfven and fast magnetosonic wave cascades are, re- 
gardless of their origin, both cut off by ion-neutral collisions 
at scales larger than the proton mean free-path, Ap, i.e., in the 
collisional regime - and thus at scales considerably larger than 
the Larmor radii of interstellar positrons. In consequence, posi- 
trons will find no Alfven or fast magnetosonic wav es from an 
MHD cascade to resonantly interact with (see also Higd on et al.l 
(l2009l) V 

The situation is completely different in the ionized phases of 
the ISM. There, the Alfven wave cascade develops with insigni- 
ficant (collisional) damping down to Ap. It then enters the col- 
lisionless range, where it is eventually cut off by linear Landau 
damping around the proton inertial length. Thus, the extended 
inertial range of the Alfven wave cascade leaves some room for 
possible resonant interactions with positrons. The fast magneto- 
sonic wave cascade, for its part, suffers strong collisional (vis- 
cous) damping. In the warm ionized phase, this damping is suf- 
ficient to kill the cascade (with the possible exception of quasi- 
parallel waves) before it enters the collisionless range. In the hot 
ionized phase, the cascade manages to reach the collisionless 
range, but it is then quickly destroyed by linear Landau dam- 
ping (again with the possible exception of quasi-parallel waves). 
Altogether, no fast magnetosonic waves from an MHD cascade 
have sufficiently small scales to come into resonant interactions 
with positrons. 



2.3. Positron interactions witli MHD wave cascades 

Our previous discussion indicates that the MHD cascades 
are truncated at scales several orders of magnitude larger than 
the Larmor radii of positrons produced by radioactive decay 
or as cosmic rays, except in the hot and warm ionized phases 
of the ISM. In these phases, Alfven wave turbulence is expec- 
ted to cascade nearly undamped down to scales close to the 
Larmor radius of MeV positrons. However, as we now argue, 
this does not necessarily mean that short-wavelength Alfven 
waves will resonantly interact with MeV positrons. Indeed, 
magnetic fluctuations are probably highly anisotropic at small 
scales, in the sense that turbulent eddies are strongly elonga- 
ted along t he mean magnetic field, or , in mathematical terms , 
k^ » /ty (iGoldreich & Sridhaj [19951: lYan & Laz^ianl |2004 . 
Because of the important anisotropy of magnetic fluctuations, 
which increases towards smaller scales, scattering off Alfven 
waves appears questionable. The elongated irregularities asso- 
ciated with anisotropic turbulence average out over a particle 
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gyration (IChandranll2000h . If some Alfven waves are present at 
scales ~ r'^, the scattering frequency is reduced by more than 
20 orders of magnitude com pared to the situation w ith isotropic 
turbulence (see, for instance. I Yan & Lazarianll2002|). Slab turbu - 
lence gives the same order of estimates ( Yan & Lazarianll2b04l) . 
In consequence, scattering off Alfven wave turbulence should 
be extremely inefficient at confining positrons in the ionized 
phases of the ISM. In this view, the diffusion model s for positron 
transp o rt adopted iii a series of recent papers by iParizot et"aLl 
(I2005h : ICheng et al.l (|2006|) : iJean et alj (|2006|) . using diff^usion 
coefficients derived from quasi-linear theory, overestimate the 
confinement by plasma waves. 

To circumvent this problem, lYan & LazarianI (|2004|) recon- 
sidered scattering off fast magnetosonic waves, emphasizing the 
isotropy of the fast wave cascade. As we saw earlier, in the hot 
and warm ionized phases, fast waves decay away at scales much 
larger than the Larmor radii of positrons from radioactivity, ex- 
cept possibly at quasi-parallel propagation. If propagation angles 
are not or only weakly randomized by wave-wave interactions 
or by chaotic divergence of magnetic field lines, quasi-parallel 
waves may survive down to much smaller scales. In that case, 
they may be involved either in gyroresonance or in Cherenkov 
resonance, also known as transit-time damping (TTD) reso- 
nance, with positrons. As explained at the end of section 12.11 
Cherenkov resonance occurs at propagation angles 0t such that 
cos ©T ^ Vms/(y cos a), independent of the wavelength. But we 
know that only those fast waves with ©x - have a chance to es- 
cape heavy damping. From this, we conclude that only positrons 
with pitch angles satisfying cosff ^ Vms/f^ 1 have a chance 
to experience TTD resonance. At precisely cosff - Vm^/u, the 
TTD mechanism vani shes, but its rate rises rapidly as cos a in- 
creases above V^s/v (ISchhckeiser & Milleilll998l) . Ultimately, 
positrons produced by radioactive decay are unlikely to effi- 
ciently interact with MHD waves from a direct cascade gene- 
rated at large scales. Such positron-wave interactions appear 
to be completely ruled out in the neutral phases of the ISM. 
In the ionized phases, they could potentially take place, but 
only under very restrictive conditions, including quasi-parallel 
fast waves (©t - 0) and nearly perpendicular positron motion 
(a ^ arccos(y,ns/u)). 

Of course, an additional local injection of MHD waves at much 
smaller scales could participate in the confinement of positrons. 
Some aspects of this possibility will be discussed in section [ 



2.4. Wave injection tiirougli plasma instabilities 

As we saw in the previous subsection, MHD waves injected 
at large scales into an MHD cascade are generally unable to effi- 
ciently interact with positrons from radioactive decay. However, 
MHD waves can be injected into the ISM by a variety of fluid or 
kinetic instabilities, which involve changes over the whole or a 
fraction of the velocity distribution of some particle population. 
These waves can be injected at scales « Linj, possibly directly 
into the coUisionless regime. Each case requires a dedicated in- 
vestigation of the wave damping process. The scale of the tur- 
bulence injection is a key parameter controlling the interaction 
between MHD waves and positrons. First, Linj controls the ani- 
sotropy of the Alfven wave cascade at the scale of the Larmor 
radius of positrons. Second, it also controls the cutoff wavenum- 
ber and the propagation angle of the fast magnetosonic cascade 
(see Eq. ( IA.18I )). In this subsection, we focus on some particular 
aspects of one type of kinetic instability. 

One of the most widely studied kinetic instabilities is trig- 
gered by the streaming of cosmic rays in the ISM, with a 



bulk velocity larger than a fe w times the local Alfven speed 
(IWentzellll974l: ISkilling|ll975h . The streaming instabihty is ex- 
pected to develop mainly in the intercloud medium. Cosmic-ray 
streaming compensates for the sink in the low-energy cosmic- 
ray population due to strong ionization losses inside mole- 
cular clouds. Low-energy cosmic rays scatter off their self- 
generat ed wav es, and are, therefore, excluded fro m molecular 
cloud s dLerche 1967; Skilling & Strong 1976; Cesarskv & Volkl 
119781; iDogel & Sharovl Il985h . This scen ario wa s adap ted to 
the transport of cosmic-ray electrons by iMorfilll (|1982|) . The 
streaming instability and other kinds of kinetic instabilities re- 
cently received new attention in the context of cosmic-ray diffu- 
sion in anisotropic MHD turb ulence ( [Farmer & Goldreichl2004l ; 
iLazarian & Beresnvakll20()6l) . 

Here, we restrict our discussion to the sole streaming instabi- 
lity. The waves are generated at scales k^^^ close to the gyroradii 
of low-energy cosmic rays, so that k^^ <g: Linj. If the waves ge- 
nerated by cosmic rays are to serve as scattering agents for low- 
energy positrons, then the Landau-synchrotron resonance condi- 
tion has to be fulfilled by both species, i.e., by virtue of Eq. ([3]) : 



Pv 



cos ffn 



COS ffel 



(19) 



Eq.(fT9]l implies that the wave-generating cosmic rays and the 
scattered positrons must have comparable momenta, unless the 
ratio of angular factors is very different from unity. Now, the 
wave-generating cosmic rays have typical momenta in the range 
[Pmin, Pmaxl- Their maximum momentum is set by the condition 
that waves can indeed be amplified by cosmic-ray streaming, 
which leads to p^ax ~ 100 MeV/c ( Skilling & Strong 197^. 
Their minimum momentum is set by the highest possible fre- 
quency of MHD waves, exactly as for positrons, and is there- 
fore given by p,^in - Va (see Eq. (|7]i). This means that, in 
the warm ionized phase, positrons with kinetic energies below 
- 2.5 - 10 keV (depending on the exact ion density) can no lon- 
ger resonate with Alfven waves generated by cosmic-ray strea- 
ming (see Table [TJ. We conclude that this mechanism does not 
operate over a broad range of positron momenta and that it pro- 
bably does little to confine positrons within the ionized phases 
of t he ISM. 

iHigdon et al.l ( l2009l) proposed an alternative mechanism 
whereby positrons scatter off their own self-generated waves. 
The growth rate of such an instability is proportional to the den- 
sity of resonant particles. The positron den sity depends strongly 
on their position with respect to the sources. iHigdon etal.l (120091) 
argued that whistler waves are heavily damped in the neutral 
phases of the ISM, but the streaming instability may also ope- 
rate close to the sources (for instance close to supernova rem- 
nants ; see Ptuskin et al. 2008) or above the Galactic disk where 
the ionization fraction can be close to unity. A complete estima- 
tion of this process deserves a detailed investigation and will be 
considered in a forthcoming paper 



2.5. Positron transport in dissipated turbulence and 
re-acceleration 

In the solar wind, MeV positrons can also interact with per- 
turbations that fall into the dissipative range of the turbulence, 
above the steepening observed at a fraction of the proton gyro- 
frequency at wavenumbers k > ^a- It has been proposed that 
low-frequency non-resonant magnetosonic waves can dominate 
the propagation of sub-Me V particles if magnetosonic w a ves are 
present in the solar wind dToptv gini 1 1 9851: iRagol l2006l) . iRagotI 
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d 19991) suggested that non-resonant fast magnetosonic waves 
can produce efficient angular scattering through pitch angles 
a ^ 7r/2 and thus govern the particle mean free -path at ener- 
gies ~ 1 MeV. This result could be applied to the stellospheres 
of massive stars (about a few parsecs) and to larger regions of 
the ISM (see discussion in lHigdon et al. I (l2009h . though in a dif- 
ferent perspective). 

The electric component of the waves induces a variation in 
particle energy. The average effect of such interactions leads to a 
stochastic energy gain, also known as second-order Fermi acce- 
leration. Particle stochastic acceleration (also referred to as re- 
acceleration) is important if the re-acceleration time is shorter 
than the energy loss time, i.e., if fstoch ^ floss- The re-acceleration 
time is related to the momentum diffusion coefficient Dp through 
fstoch = IDp (here, we consider only pitch-angle averaged 
quantities). It can also be expressed in terms of the particle angu- 
lar diffusion frequency Vj as fstoch - (c/Va)^ ^s^'- Then the above 
timescale ordering is equivalent to 



(20) 



Particles with energies < 1 GeV mostly suffer ionization and 
Coulomb losses. The energy loss rate depend s on whether the 
medium is ionized or neutral. lGinzburg|(ll979l) gives 



1 

Aoss 

in an ionized medium and 



73.6 -Hln(y/We.cm-3) 



1 

^loss 



20.5 + ln[(r-l)(r--l)l 



y6(r-l) 



(21) 



(22) 



in a neutral medium, where crj is the Thomson cross-section, 
jS = vjc and, as before, y is the particle Lorentz factor, «e the 
free-electron density and hh" the density of hydrogen atoms. The 
condition ( |20l ) translates into an upper limit on the mean free- 
path ; 



A < (3.4 X 10" V) 



ln(y/«e.cm-3) + 73.6 



(23) 



3. Effects of collisions with gas particles 

In this section, we examine the case when the trajectories of 
positrons in the ISM is only driven by their collisions with gas 
particles while they propagate along a steady state magnetic field 
line, without considering interactions with MHD waves, an as- 
pect usually overlooked by the previous analysis. During these 
collisions, high-energy positrons not only lose energy, but they 
also undergo pitch angle scattering. The scattering of kinematic 
parameters of positrons reduces the maximum distance they can 
travel along a straight line in the ISM. In section [TT] we review 
and describe the interaction processes between positrons and in- 
terstellar matter In section [T2l we describe the methods used to 
calculate positron propagation in this so-called "collisional re- 
gime", and in section [331 we present the detailed results of our 
computations for 1 MeV positrons. The main results obtained 
for positrons with initial kinetic energies ranging from 1 keV to 
10 MeV in typical ISM phases and the case of propagation in a 
turbulent magnetic field are described in section [ 



in an ionized medium and 



3.1. Positron interactions with interstellar matter 

In the absence of collisions, positrons move along magnetic 
field lines in helical trajectories. When positrons interact with 
gas particles, they can either gain or lose energy in elastic and 
inelastic scattering or even annihilate with free or bound elec- 
trons. The energy and pitch angle variations resulting from the 
interaction depend on the energy of the incident positron and on 
the velocity and nature of the target particle. 

The positrons that we are studying spend most of their life- 
time travelling at high energy, with a kinetic energy greater than 
the thermal energy of target particles in the ISM. Consequently, 
we do not account for the propagation of positrons at thermal 
energy, but postpone the discussion of this case to section [3T4l 

Interactions of positrons with the ionized component of the 
interstellar gas and, at high energy (E > 100 MeV), with the 
magnetic field (synchrotron radiation) and the interstellar ra- 
diation field (inverse Compton scattering) are generally consi- 
dered as continuous processes. They are referred to as conti- 
nuous energy-loss processes and are presented in section [3.1.11 
In a mostly neutral gas, positrons lose a larger fraction of their 
energy when they excite or ionize atoms or molecules. Such in- 
teractions, which result in a quantified variation of the positron 
energy, are described in section [J.1.2| 



^<(3.4x lO^V) 



«H» cm-3 ni,( 



(r-DyS' 



In [(r - l)(r' - D] + 20.5 
(24) 

in a neutral medium. This implies that if, in a particular me- 
dium, the angular diffusion frequency is high enough, particles 
are reaccelerated. In that case, a solution can be obtained with 
the help of a diffusion-convection equation. This type of investi- 
gation is beyond the scope of the current paper and is postponed 
to a future work. However, the above discussion a lready calls 
some of the results derived by iHigdon et al.l ( l2009l) into ques- 
tion, as these authors ignored re-acceleration in their analysis. 
What emerges from the above discussion is that re-acceleration 
can be important, especially in low-density ionized phases. More 
specifically, if we i ntroduce the parameter values proposed by 
iHigdon et al.l ( l2009l) for the Galactic interstellar bulge (see their 
Table 1), we find that positron re-acceleration is important in the 
very hot phase of the inner bulge, and probably in the hot phase 
of the middle and outer bulge. 



3.1 .1 . Continuous energy-loss processes 

Figure[T]shows the energy loss rates of positrons as functions 
of their energy in a fully ionized plasma with temperature T = 
8000 K. The energy loss rates by synchro tron radiation and in- 
verse C ompton scattering are derived from [Blumenthal & Gouldl 
( [1970|) . They are proportional to the magnetic-field and photon 
energy densities, respectively. The energy loss rate by brem- 
sstrahlung is cal culated using the approximation presented in 
iGinzburd (119791) . The energy loss rate due to Coulomb colli- 
sions (free-f ree) is calculated u sing the Bhabha cross-section as 
described in [Asano et an (l2007h . 

In a fully ionized plasma. Coulomb collisions represent the 
dominant loss process for positrons with £ < 100 MeV. Since the 
energy loss rates by synchrotron radiation and inverse Compton 
scattering depend on the magnetic-field and photon energy den- 
sities, they vary with the location of positrons in our Galaxy. It is 
expected that inverse Compton losses become dominant in zones 
close to stellar clusters. However, we are interested in positrons 
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Fig. 1. Energy-loss rates of positrons in a 8000 K plasma with 
number densities ris - «h+ - 1 cm"-'. The long-dashed line 
shows the energy-loss rate due to synchrotron radiation assu- 
ming Bo - 5 ijG and an average value of sin^ff = 2/3. The 
dotted line is the energy-loss rate by inverse Compton scattering 
assuming Uph - 0.26 eV cm"-', which corresponds to the cos- 
mic microwave background energy density. The dot-dashed line 
shows the energy-loss rate by bremsstrahlung emission due to 
collisions with charged particles. The solid line corresponds to 
the positron energy-loss rate due to Coulomb interactions bet- 
ween positrons and electrons. 



with kinetic energy < 10 MeV and, therefore, we may neglect 
synchrotron losses as long as B < 1 mG. 

The average scattering angle induced by Coulomb collisions 
occurring over a time interval St is estimated through the relation 



sin 6 •■ 



t/<shr0) 
^St 



dt 



(25) 



where d{sin^ 6)/dt is the rate of variation of the average sin^ 6 
and 6 is the positron scattering angle in the laboratory frame. 
The rates of variation of the average scattering angle are evalua- 
ted by integrating the Bhabha (e+e" collision) and Rutherford 
(e^p collision) d ifferential cross sections (see Appendix in 
lAsano et al.ll2007l) . over the scattering angles as 



d{sm 



dt 



f 



sin^ e ■ 



do-* 



dcos 0* 



sm6*de* 



(26) 



where n is the number density of target particles, 6* is the scatte- 
ring angle in the center-of-mass frame, „^ is the differential 
cross section of the intera ction and 0* is the minimum scatte- 
ring ang le (see Eq. (B7) in lDermeiill985h . 

While the energy loss rate from e^p collisions is negligible 
compared to that from e^e" collisions, their deviation rates are 
equivalent. We checked that the resulting deviation rates are in 
ag reement with those derived at low energy using the formalism 
of lHubal(l2006h . 

3.1.2. Anniliilation and interactions witli atoms and 
molecules 

In a neutral medium, high-energy positrons lose energy 
mainly by ionizing and exciting atoms and molecules, or 
they annihilate directly with bound electrons, the elastic 



scattering process being ne ghg ible at these ener gies (see 
ICharlton & HumberstonllfoOOl and lWallvn et al.lll994l) . Such in- 
teractions occur at random while positrons propagate in the ISM. 
Therefore, we will determine the kind of interaction and cal- 
culate the variation of the kinematic parameters with a Monte- 
Carlo method that incorporates the corresponding cross sections. 

The ionization and excitation cross sections as well as the 
differential cross sections as functions of the ener gy lost by 

J ositron s in ionizing collisions were calculated by iGrvziriskil 
1965allbll3) . The cross section of annihilation in flight of posi- 
trons with bound electrons is equal to that of annihilation with 
free electrons, since the binding energy of electrons is negligible 
with respect to the kinetic energy of the positrons under study. 
We use t he cross section of anni hilation with free electrons pre- 
sented in lGuessoum et al.l (l2005l and references therein). 

In the Monte-Carlo simulations, the energy lost by a positron 
when it ionizes an atom/molecule is chosen randomly according 
to its differential cross section. The energy lost by a positron 
when it excites an atom/molecule is derived from the energies 
of the atomic levels involved in the interaction. Once the energy 
lost by ionization or excitation is known, the scattering angle of 
the positron is calculated with the kinematics of the interaction, 
assuming that the atoms/molecules stay at rest and assuming azi- 
muthal symmetry. 

3.2. Simulations of the collisional transport 

Our Mo nte-Carlo simulation s are based on the method s 
presented in iBussard et al.l (Il979h and iGuessoum et all (l2005h . 
However, we add the calculation of the trajectories of positrons 
while they propagate along magnetic field lines and include ad- 
ditional steps to account for the scattering of positrons by colli- 
sions. 

At the initial time {k = 0), positrons are located at (x, y, z, t) = 
(0, 0, 0, 0) in a frame such that the magnetic field Bq is directed 
along the z axis. The initial kinetic energy of positrons is E^. The 
direction of their initial velocity i»o, which is defined by their 
initial pitch angle ao and their initial phase Oq (see Figure |2]i, is 
chosen randomly according to an isotropic velocity distribution 
in the entire space. 

We proceed through successive iterations. At step k, the re- 
sulting pitch angle (ait) and phase (O^) are given by 



cos at - cos 1 COS0 + sino-i:-! sin cost 



COS((Di - <t>A._i) 



COS0- COSQ;;t_l COSQ-J: 



(27) 



(28) 



smort-i sino'^. 

where the azimuthal scattering angle ^ is a random number uni- 
formly distributed between and 27r. The kinetic energy Ek after 
collision w ith an atom/molecule or annihilation with an electron 
IS given bv lBussardetaLl(ll979l) : 



R = exp 



dEJdt 



dE 



(29) 



where R is a random number uniformly distributed between 
and 1, a-j j is the cross section of process / (see section [3.1.2b 
with target j and nj is the number density of target j. The kind 
of interaction is chosen randomly according to its probability 



nj o-jjiEt) 
Z n j o-ij(Eit) 



(30) 
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The energy lost and the scattering angle induced by collisions 
are derived using the methods presented in section [J.1.2| 

The increment in time between two iterations, 6t, is equal to 
a scattering time Sts defined as 



6ts = 



d{wte)ldt 



(31) 



with (i(sin^ 6)ldt the rate of variation of the average sin^ (see 
Eq. |26] l and e a prescribed number <«; 1 . During <5fj, positrons 
lose an amount of energy dE/dtxdts via continuous energy-loss 
processes, and, they are slightly scattered by an angle 6 obtained 
from Eq. |25] If the time interval between two successive colli- 
sions with neutrals is shorter than the scattering time 6ts (e.g., 
in a weakly ionized medium), then the increment in time 6t is 
replaced by the time interval between two collisions. 

Between two steps and/or two interactions, positrons propa- 
gate in a regular helical trajectory along magnetic field lines with 
a gyroradius r'^ sin ak-i (defined in section 12.1b . The positron 
position after each step is therefore calculated analytically. 

The above procedure is repeated until positrons annihilate or 
their kinetic energy falls below 100 eV, since the distance trave- 
led below this energy becomes negligible. To estimate the spatial 
distribution of positrons and their lifetime, we read out the po- 
sitron location and slowing-down time at the end of the track. 
To save CPU time while obtaining sufficiently accurate results, 
Monte-Carlo simulations are performed with a number of posi- 
trons ranging from 5000 to 20000. 

To validate our code, we performed several tests for posi- 
trons with Eq < 10 MeV released in different media and we 
compared the results of our simulations with previous work. 
The fractions of positrons annihilating in flight with free or 
bound electrons are in agreement with the results presented by 
[Beacom & Yiiksel (2006) and Sizun et al. (2006). At low energy 
(£() < 1 keV), the fractions of positronium formed in flight by 
charge exchange with atoms and molecules in different media 
are identical to those obtained bv lGuessoum et al" I (l2005h . Once 
positrons are thermalized, their propagation behaves as classical 
diff\ision. The distance they travel then is negligible compared to 




Fig. 2. Kinematic parameters of the scattering in the laboratory 
frame. The magnetic field Bq lies along the z axis, ita-i is the 
pitch angle between the positron velocity Uk-i and Bo before the 
interaction, is the pitch angle after the interaction, when the 
positron velocity is U/.. <!>/._ i and <^i, are the phases of incident 
positrons before and after the interaction, respectively. is the 
scattering angle of positrons and (p is the azimuthal scattering 
angle. 



the distance traveled in the slowing-down regime, except in the 
hot medium (see section U4l ). 



3.3. Results of the simulations for 1 MeV positrons 

This section presents, as an illustrative example, the detailed 
results of simulations obtained for positrons with initial kinetic 
energy £0=1 MeV. For this first set of simulations, we chose 
the warm medium, whose temperature is intermediate between 
those of the cold and hot media. The warm medium is also the 
most appropriate medium to study the impact of the ionization 
fraction, which covers the whole possible range between and 
1. We assumed the magnetic field to be uniform, with a strength 
Bo - 5 fiG. 

Since the energy loss rates and the frequency of inelas- 
tic interactions are proportional to the target density, for given 
abundances and a given ionization fraction, the distance tra- 
veled by positrons scales as the inverse of the total density. 
Therefore, we performed simulations for «h - 1 cm"^, with 
«H - «H» + «H+ + 2nH, the total number density of hydro- 
gen nuclei. Preliminary tests of our simulations indicate that the 
presence of ionized or neutral helium with an abundance ratio 
He/H X! 0.1 does not strongly affect positron propagation (dif- 
ferences less than 10%). We also neglected molecular hydrogen, 
except in the MM where H2 molecules are the dominant species. 

Figure[3]shows the spatial distributions of positrons along the 
magnetic field direction (z), as extracted from our simulations, 
once they reach E - 100 eV. Positrons are initially "injected" 
at (x, y, z) = (0, 0, 0) with a pitch angle chosen randomly ac- 
cording to an isotropic velocity distribution in the entire space. 
The temperature is 8000 K and the ionization fraction ranges 
from to 1. The field-aligned distributions are nearly uniform 
out to the maximum distance traveled along field lines, c/max (ob- 
tained when the pitch angle is and remains equal to 0). This 
means that the pitch angle does not change significantly over 
most of the slowing-down period of positrons. If the pitch angle 
remained strictly constant, the distance traveled along field lines 
would be equal to dmax x cos ao. Since positrons are initially 
emitted isotropically, cosao is uniformly distributed between 
-1 and 1, so that the field-aligned distributions of positrons at 
the end of their slowing-down period would be uniform. In rea- 
lity, the slight pitch angle scattering induced by collisions pro- 
duces a slight scattering of their final field-aligned positions and, 
therefore, smoothes out the edge of the otherwise uniform dis- 
tributions. The maximum distance d^^x is obtained by integra- 
ting over energy the ratio of the positron velocity to the energy 
loss rate. In the present case, the energy loss rate is the sum 
of the contributions from Coulomb collisions (see section [J. 1.1b 
and inelastic interactions (ionization and excitation) with atoms 
and molecules. The l atter contribution s are evaluated using the 
Bethe-Bloch formula dOinzburd 19791) . 

The extent of the spatial distributions in the directions per- 
pendicular to the magnetic field is a few times the Larmor radius, 
i.e., negligible with respect to the extent of the field-aligned dis- 
tributions. We characterize the extent of the field-aligned distri- 
butions presented in Figure |3] by their full width at half maxi- 
mum (FWHM//). Its uncertainty is calculated with a bootstrap 
method taking into account the uncertainty in the number of po- 
sitrons (Poisson statistics) in each spatial bin. 

The extent of the field-aligned distributions increases natu- 
rally with decreasing ionization fraction (see Figure [3), since 
the energy losses through Coulomb collisions dominate over 
the losses due to inelastic interactions (ionization and excita- 
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tion) with neutrals. Depending on the ionization fraction, the half 
FWHM// is 15% to 30% lower than d^^^. 

3.4. Transport of positrons in the different ISM piiases 

This section summarizes the results of the simulations for 
positrons with initial kinetic energies ranging from 1 keV to 
10 MeV, in the different ISM phases. 

Figure|4]shows the extent (FWHM//) of the field-aligned dis- 
tribution of positrons at the end of their slowing-down time (i.e., 
when they reach 100 eV), as a function of their initial kinetic 
energy in the different ISM phases. As in section [33] , the hy- 
drogen density is arbitrarily set to hh = 1 cm""*, and the ioniza- 
tion fraction is 0, 0.001, 0.1, 0.9 and 1 in the MM, CM, WNM, 
WIM and HIM, respectively. The extents presented in Figure |4] 
are compared to the extents that the distributions would have if 
the positron pitch angles remained constant, namely, 2dmax- In 
each ISM phase and for any given initial energy, FWHM// is si- 
milar to, and always slightly less than, 2d„iax, as expected from 
the effect of pitch angle scattering. 

The above FWHM// extents are given assuming a uniform 
magnetic field in the z direction. However, field lines in the 
ISM are perturbed by turbulent motions. As a result, realistic 
magnetic fields in the ISM consist of a mean ("regular") com- 
ponent plus a turbulent component, leading to chaotic field lines. 
The distances traveled by positrons along the uniform field (i.e., 
along the z axis) are shorter than the total distances traveled 
along the actual chaotic field lines. 

We estimate the effects of the turbulent magnetic field 
on the distances traveled by positrons by adding a turbulent 
component to the uniform field. This turbulent field is mo- 
de lled using the plane wave a pproximation method presented 
bv iGiacalone & Jokipiil ( |1994|) . We assume a ratio 6B/B0 - 1. 
The fluctuations follow a Kolmogorov spectrum and have a 
maximum turbulent scale /Imax - 10-lOOpc in the hot and warm 
phases and /l„,ax - 1-10 pc in the cold neutral and molecular 
phases. These scale lengths correspond to the typical sizes of the 
respective phases. Since the Larmor radii of positrons with ki- 
netic energies in the considered range are extremely small com- 
pared to the turbulent scale lengths, positrons simply propagate 
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Fig. 3. Spatial distributions of positrons along field lines at the 
end of their slowing-down period in a warm medium {T - 8000 
K and «h = 1 cm "*), for several ionization fractions, /ion- The 
distributions are normalized to unity and the error bars are esti- 
mated according to Poisson (counting) statistics. 
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Fig. 4. FWHM of the field-aligned distributions of at the end of 
their slowing-down period, as functions of their initial kinetic 
energy. The hydrogen density is arbitrarily set to «h = 1 cm"^ 
in each ISM phase. The lines show the distance 2(imax (see text 
for details) while points show FWHM// of the distribution as 
obtained by the simulations. 

along the turbulent field lines. Under these conditions, the ac- 
tual coordinates of positrons at the end of their slowing-down 
period are obtained by carrying their coordinates along the uni- 
form field (as calculated in section |3.3l l over to the curvilinear 
frame of the turbulent field lines. These coordinates are calcu- 
lated by Monte-Carlo simulations, using a large number of ran- 
domly chosen turbulent configurations for a given /Imax- Figure |5] 
shows examples of the positions of positrons at the end of the 
slowing-down time, calculated with this method in a WIM with 
«H = 0.2 cm"-' and fio„ - 0.6. 

Adding such a turbulent magnetic field in the simulations re- 
duces the FWHM of the positron field-aligned distributions by 
a factor ^ 0.75 and broadens their transverse (i.e., perpendicu- 
lar to Bo) distributions due to the chaotic behavior of field lines. 
At low energy, when the distance traveled by positrons is shorter 
than the maximum scale of the turbulence, /Imax (e-g-, Figure|5^), 
positrons have a quasi-uniform distribution along any field line 
that can be considered, at this scale, as straight and randomly til- 
ted with respect to Bq. In this case, the spatial distribution of po- 
sitrons is nearly spherically symmetric. At intermediate energy 
(e.g.. Figure |5})), positrons have a more collimated distribution 
with a transverse dispersion that increases with distance from 
the source. At high energy, when the distance traveled by po- 
sitrons is large compared to /Imax (e-g-. Figure |5};), the positron 
distribution is highly collimated along Bo- 

The shape of the transverse distributiorQ is strongly peaked, 
as illustrated in Figure |6] for positrons with Eq - 1 MeV pro- 
pagating in a WIM. The transverse distribution is close to a de- 
caying exponential function. Its extent increases with /Imax, as 
expected since a larger scale length leads to larger transverse ex- 
cursions of field lines away from the uniform field. In contrast, 
the extent of the field-aligned distribution does not change with 
'imax, because, on average, positrons propagate along the uniform 
magnetic field. 

Due to the cusped shape of the transverse distribution, 
FWHMj^ is not necessarily a very meaningful quantity, as the 
fraction of positrons inside FWHM^ is low. Moreover, this frac- 

3 The transverse distribution is ^ with dS = Inrdr and r the radius 
in the xy plane. 
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tion changes significantly with A^iix (e.g., from ^1% to ^9% for 
the cases presented in Figure|6]l and with £0 (e.g., from ^1% to 
^30% for £0 decreasing from 10 MeV to 1 keV). In contrast, 
the fraction of positrons inside the FWHM of the field-aligned 
distribution ranges from ~ 80% to ~ 95%. 

Consequently, instead of characterizing the spatial extents of 
the positron distribution by its FWHM, we will use the field- 
aligned length (2z9o) and the diameter (Iryo) that contain 90% 
of the positrons. Figure |7] shows the extents of the spatial dis- 
tribution as functions of the initial kinetic energy of positrons. 
The physical conditions (density, ionization fraction, tempera- 




ture,. . . ) of the different ISM phases correspond to the extreme 
values presented in Table[T] The minimum and maximum spatial 
extents are obtained with the maximum and minimum densities 
and ionization fractions listed in Table [T] and with the minimum 
and maximum /imax, respectively. 

The ratio rgo/zgo increases up to ~0.8 with decreasing energy, 
confirming that the positron distribution is more isotropic at low 
energy. A spherically symmetric distribution would have produ- 
ced a ratio rgo/zgo = 1.2. 

The results displayed in Figure [T] show that MeV positrons 
travel long distances before reaching 100 eV. Positrons with 
£0 ^ 30 keV generally cannot escape molecular clouds (assu- 
ming the smallest size of molecular clouds is ~1 pc). On the 
other hand, positrons with Eq >1 keV generally escape the hot 
medium (assuming the typical size of hot regions is ~200 pc). 

So far, we have discussed the distance traveled by positrons 
over their slowing-down period, which ends when their kinetic 
energy falls below 100 eV. We now consider the distance they 
travel before annihilating with electrons. In the neutral phases, 
low-energy positrons {E <100 eV) annihilate mostly by forming 
positroniums through charge exchange processes with atoms or 
molecules. The distance they travel before annihilating is ne- 
gligible compared to the distance traveled at high energy (see 
Table 4 of Jean et al. 2006). The situation is different in the io- 
nized phases. There, positrons quickly thermalize and then anni- 
hilate with electrons (via radiative recombination or direct anni- 
hilation). Since the time scale of annihilation (rann) is inversely 
proportional to the electron density, positrons in low-density me- 
dia may travel a long distance before annihilating. The r.m.s. 
distance traveled by thermalized positrons along a uniform ma- 
gnetic field is 

Zth = y/2D\\T^nn , (32) 

with the diffusion coefficient - AsVe/3. Here, is the posi- 
tron r.m.s. velocity, A^^ is the positron collisional mean free-path, 

T- 

Ae = (2.5 X 10' Vm) — , (33) 




Fig. 5. Positions of positrons at the end of their slowing-down 
time in a WIM taking into account colUsional transport in a tur- 
bulent magnetic field {SB/Bq = 1 and /l^jax = 10 pc). The uni- 
form magnetic field lies along the z axis. Positrons are injected at 
the origin of the frame with initial kinetic energies of (a) 5 keV, 
(b) 100 keV and (c) 5MeV. 
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Fig. 6. Transverse distributions of positrons at the end of their 
slowing-down period in a WIM, calculated for maximum turbu- 
lent scales /Imax - 10, 100 and 1000 pc. The initial kinetic energy 
of positrons is Eg = 1 MeV. The distributions are normalized to 
unity. 
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and va lues of Tann are derived from Table 3 of iGuessoum et all 
(l2005h . 

The r.m.s. distance Zth characterizes the width of the 
Gaussian field-aligned distribution of thermalized positrons at 
the time of annihilation, when they are initially injected isotro- 
pically at z = 0. Taking into account the turbulent magnetic field, 
the 3-D spatial distribution of the sites where positrons annihi- 
late, is calculated by carrying their coordinates along the uniform 
field over to the curvilinear frame of the turbulent field lines. We 
then derive the field-aligned and transverse extents of this dis- 
tribution (2zth,90 and 2r±()o, respectively), defined as the field- 
aligned length and the diameter that contain 90% of the annihi- 
lating positrons. In the WIM, Zth90 <0.015 pc, i.e., <K Z9q- In the 
HIM, z,h,90 ~ 300 and -600 pc,' for «e = 0.01 and 0.005 cm-\ 
respectively. These values are > zgo for positrons with £0 5 
10 keV (see. Figure|7]l. Furthermore, Zth,90 in the HIM is larger 
than the typical s ize of hot regions, confirming the estimates of 
iJean et all (|2006|) that positrons do not annihilate in the HIM. 

Combining these results with the half extents (zyo and rgo) 
presented in Figure|2l we calculate the actual full half-extents of 
the spatial distribution of annihilating positrons, as functions of 
their initial kinetic energy : 
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Fig. 8. Minimum and maximum extents of the spatial distribu- 
tions of annihilating positrons, along (top) and perpendicular 
(bottom) to the uniform magnetic field, taking into account the 
turbulent behavior of the field lines as well as realistic values for 
the density in each ISM phase. 
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Fig. 7. Minimum and maximum extents of the spatial distribu- 
tions of positrons reaching 100 eV, along (top) and perpendicu- 
lar (bottom) to the uniform magnetic field, taking into account 
the turbulent behavior of the field lines as well as realistic values 
for the density in each ISM phase. 



Zann,90 = ^Zgo + ^1,90 (^4) 

and 



'•ann,90 = y'-go + ''th.go " ^^^^ 

Figure[8]presents these extents in the different phases of the ISM. 
They are similar to the results presented in Figure|7l except in the 
HIM where the lifetime of thermalized positrons is so long that 
they have time to travel a non-negligible distance before they 
annihilate. 

4. Advection with large-scale fluid motions 

As discussed in section [31 low-energy positrons are at least 
loosely coupled to the thermal plasma, regardless of their pro- 
duction channel (/?^ decay or cosmic-ray interactions). More 
specifically, positrons are tied to magnetic field lines and the lat- 
ter are frozen-in into the plasma. If the plasma is subject to orde- 
red or stochastic large-scale motions, positrons will be dragged 
along with the plasma and the field lines. 

The Galactic center region is the site of intense massive star 
formation, as recently confirmed by long-exposure Chandra ob- 
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servations dMuno et al.ll2006b ; it harbours within the innermost 
30 pc three bright and young star clusters : the Arches, the 
Quintuplet and the central Galactic cluster close to Sagittarius 
A*. Most massive stars are grouped into clusters and/or as- 
sociations of up to several hundreds. The powerful winds 
from these stars together with their terminal supernova explo- 
sions act colle ctively to produc e hot, low-densit y superbubbles 
dMcCrav & Ka fatos 1987; Par izot et al.l l2004bl) . These super- 
bubbles have a profound impact on the structure of the ambient 
ISM. Some of them are even powerful enough to break through 
the Galactic disk and produce chimneys, thereby ve nting their 
energy into the halo (iMcClure-Griffiths et al.l 120061 and refe- 
rences therein). 

Another possibility to drive strong winds in the Galaxy 
is again related to the streaming instability triggered by cos- 
mic rays escaping the Galactic disk into the halo. This phe- 
nomenon may be rather effective at distances ~ 1 kpc from 
the Galactic cent er, at the transition between the diffusion and 
advection zones jPtuskin et al.lll997l) . In reality, the transition 
could well be closer to the Galactic center. Recent HESS obser- 
vations tend to support the idea that the cosmic -ray energy den- 
sity around the Galacti c center is higher than in the solar system 
dAharonian et al.ll2006l) . 

The two next subsections examine some of these aspects in 
more detail. 

4. 1 . Massive-star forming regions : small-scale advection 

Positrons produced by radioactive decay in massive stars are 
typically injected into the hot interior cavity of a superbubble. In 
this hot medium, most positron transport takes place during the 
slowing-down period, a phase during which positrons are subject 
to collisions with gas particles. 

In addition, positrons are advected with the motions of 
the plasma filling the superbubble. Because of the presence of 
strong stellar winds and supernova shocks inside the superbubble 
and because of their interactions with each other as well as 
with local density inhomogeneities, the hot cavity is undoub- 
tedly highly turbulent (Bvkov 2001, and references therein). 
According to iBvkovl dlOOlh 's model of inhomogeneous turbu- 
lence in superbubbles, positrons, which follow field lines, are 
subject to the same chaotic transport, characterized by a diffu- 
sion coefficient D = Ui-^^ L. With a turbulent velocity in the 
range v^ms ~ (100 - 1000) km s"' and a mean scale of turbu- 
lent eddies in the range L ~ (1 - 10) pc, corresponding to the 
typical distance between two stars (and two shocks), we find 
D ~ (3 X 10-^ - 3 X 10") cm- s"'. It then follows that over 
their lifetime (rnfe ~ (10^ yr)/«H,cm-^ with «h given in Table [1]; 
see sectionO, positrons diffuse over distances ~ (80 - 1000) pc. 
Thus chaotic motions appear to "confine" low-energy positrons 
in the hot medium rather efficiently compared to collisions and 
anisotropic MHD turbulence. 

Large-scale MHD compressible fluctuations also provide 
stochastic acceleration. The rate of re-acceleration is given by 
Tacc ~ i^rms/9i in the limit u^^sL > vA^fp, where v is the 
positron velocity and A^fj, its mean free-path controlled by 
pitch-angle scattering (Bykov 2008, private communication). 
Re-acceleration is efficient only if the re-acceleration rate is large 
compared to the Coulomb energy loss rate (see Eq. |23] ). i.e., if 
I'rms > (80 kms~') (L/10 pc) «e,cm-' y6 ' ■ Re-acceleration can 
contribute to delaying the annihilation of positrons. 

At the scales explored by the gyromotion of positrons pro- 
duced by radioactive decay or by cosmic-ray interactions, one 
should also account for transport processes similar to those in- 



voked for low-energy particles in the solar wind. iRagotI ( Il999l) . 
for instance, considered the role of non-resonant fast magneto- 
sonic waves in the dissipation regime of turbulence, i.e., above 
the steepening observed at a fraction of the thermal proton gyro- 
frequency (see lAlexandrova et al.ll2008l and references therein). 
This non-resonant transport process can produce efficient angu- 
lar scattering through pitch angles a n/2 and thus govern the 
particle mean free-path at energies ~ 1 MeV. We believe that 
this process adapted to the stellospheres of massive stars (typi- 
cally a few parsecs in size) should be relevant for the problem of 
positron transport. 

A full investigation of the local transport of positrons within 
the stellospheres of massive stars and of their spatial transport 
and re-acceleration in chaotic compressible turbulence is beyond 
the scope of this paper and will deserve future separate studies. 

4.2. Galactic winds : large-scale advection 

Galactic winds can produce a systematic shift of the fluid at 
velocity y„ toward high Galactic latitudes. Positrons coupled to 
the fluid will be transported over their lifetimes Tjife over typical 
distances ~ y„ Tiife (for a constant wind velocity). As mentioned 
above, large-scale winds can be driven termally by the collective 
action of massive stars (stellar winds and supernova explosions) 
in active regions of the Galaxy or by the streaming of cosmic 
rays into the Galactic halo. 

Several observations support the idea of starburst episodes 
in the histor y of our Gala xy (some of their properties were 
reviewed bv IVeilleux et al.l 12005). There are further pieces of 
evidence indicating an active nuclear wind in our Galaxy. 
For instance, powerful mass ejections from the Galactic cen- 
ter were observed on scales from a few arcminutes to tens 
of degrees at several w avelengths (iMorris & Serabvnl 119961 : 
lYusef-Zadeh et al.ll2000l) . Some evidence for a large-scale bi- 
polar win d from the Galactic center was f ound in infrared dust 
emission ( iBland-Hawthorn & Cohenl2003h . The authors attribu- 
ted their 8.3 fim observations to dust entrained in a Galactic wind 
that was powered by a central starburst several million years 
ago. Wind velocities in the Galactic ce nter are difficult to esti- 
mate. iBland-Hawfliorn & Cohen! (l2003b considered vertical ve- 
locities ~ 150 km s ' at Galactic hei ghts of a few 100 p c. A 
recent HST and FUSE data analysis ( iKeenev et al.l 120061) un- 
covered four high-velocity absorption components with typical 
velocities ~ 250 kms"' at a distance of a few kpc from the 
Galactic plane. One of these components was blueshifted and 
interpreted as being part of a Galactic fountain returning mat- 
ter toward the Galactic center. Alternatively, these observations 
might correspond to wind velocities at high altitudes (> 1 kpc) 
that would be driven by the combined effect of the flux of cos- 
mic ray s escaping into the halo and the pressure of the thermal 
plasma (lEverett et al.ll2008l) . 

The wind closer to the Galactic disk or the Galactic c enter is 
certainly slower, and it probably speeds up with height. iTotanil 
(l2006b . instead of a starburst episode, favored a wind with ve- 
locity ~ 100 km s ' produced by a flaring activity from Sgr A*. 
Assuming convective transport (no other type of transport was 
considered in his analysis) throughout the lifetimes of positrons 
makes it possible to fill the extended 51 1 keV bulge. 

Another constraint on wind velocities, probably more rele- 
vant to positrons produced in the Galactic disk, comes from the 
presence of cosmic-ray radioactive elements and the analysis of 
secondary-to-primary ratios. Incorporated into a transport mo- 
del, these observables are best fitted with a wind of constan t 
velocity < 15 km s ' (IStrong et al.l 120071: iMaurin et al.ll2002l) . 
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Under the assumption of a constant wind velocity of 15 km s ' 
and for typical positron lifetimes ~ (10^ yr)/«H,cm-3 (see sec- 
tion |3]l, the distances covered by convection only are rather mo- 
dest (~ (1.5 pc)/«H,cm-3)- We do not expect this effect to be do- 
minant in the transport of positrons from the Galactic disk. 



5. Conclusion 

The purpose of this work was to carefully examine the pro- 
pagation of positrons with kinetic energies <10 MeV in the dif- 
ferent phases of the ISM. We identified and analyzed three main 
transport mechanisms : scattering off MHD waves, ballistic mo- 
tion along magnetic field lines and advection with large-scale 
fluid motions. 

In section [2] we presented an investigation of the scattering 
off" of positrons by MHD waves. Two necessary conditions are 
required for this process to work : (1) MHD waves and positron 
motions must satisfy a synchrotron resonance condition and (2) 
the corresponding MHD waves must exist. Since MHD waves 
exist only at frequencies lower than the local proton cyclotron 
frequency, the resonance condition and, therefore, the scattering 
off" MHD waves are restricted to positrons with kinetic energies 
higher than a threshold which depends on the physical parame- 
ters of the local ISM (see Table We showed that positrons 
in the neutral phases (i.e., WNM, CM and MM) do not scatter 
off MHD waves, because the latter are damped by collisional ef- 
fects, causing MHD cascades to be truncated at scales several 
orders of magnitude above the Larmor radii of MeV positrons. 
In the ionized phases. Landau damping cuts off the Alfven cas- 
cade at the dissipation scale. In principle, positrons could be in 
resonance with this part of the MHD cascade. However, the ani- 
sotropy of magnetic perturbations in the ISM leaves only a small 
fraction of the turbulent energy injected at large scales, thereby 
casting some doubt on this conclusion. 

We also briefly discussed the possible action of resonant 
Alfven waves generated by the streaming of cosmic rays in the 
WIM. It appears that positrons can resonate with these waves 
over a rather restricted momentum range. In the case of self- 
generated streaming modes, the issue is uncertain : while the 
waves are probably rapidly damped in a neutral medium, the 
interaction of positrons with the waves injected in an ionized 
medium may well help to confine them near their sources over 
longer timescales. 

We also addressed the effect of compressional magnetosonic 
waves produced either in the large-scale cascade or injected lo- 
cally, as in the solar wind. These waves can indeed produce some 
particle transport and re-acceleration, especially in the HIM. The 
importance of re-acceleration requires a kinetic treatement of the 
transport problem,i.e., solving a diffusion-convection equation. 

Positrons that do not scatter off MHD waves move along ma- 
gnetic field lines with helical trajectories which are perturbed by 
collisions with the neutral and/or charged particles of the ISM. 
A detailed analysis of this propagation mode, also called "colli- 
sional regime", is presented in section[3] We found that the pitch 
angle of positrons is only slightly perturbed by collisions during 
their propagation at high energy. The spatial distribution of posi- 
trons at the end of their lifetimes was calculated using a Monte- 
Carlo method, assuming that positrons are isotropically emit- 
ted at a point source in the different ISM phases, with various 
characteristics of the turbulent magnetic field. The spatial dis- 
tributions of high-energy positrons are nearly uniform along the 
mean magnetic field, but with an extent ~1.7 times smaller than 
twice the maximum distance traveled by positrons. Turbulence 



induces a scattering of positrons perpendicular to the mean ma- 
gnetic field, which cause an extent of the transverse distributions 
that is generally negligible compared to the extent of the field- 
aligned distributions, except for positrons at low energy (E<10 
keV) or in high-density media (e.g., in the MM). The distances 
traveled by thermal positrons before they annihilate are negli- 
gible compared to the distances traveled at high energy, except 
in the HIM. There, their lifetimes are sufficiently long to allow 
them to diffuse along field lines over distances larger than the 
typical size of the HIM. 

In this collisional regime, positrons can travel over long dis- 
tances (e.g., < 30 kpc/nn.cm ' for £0 = 1 MeV). Consequently, 
they go through several phases of the ISM before annihilating. 
In this case, since positrons probably propagate along field lines, 
the spatial distribution of their annihilation sites should depend 
on the spatial distribution of the Galactic magnetic field. 

Finally, in section [H we reviewed the process of transport 
by large-scale fluid motions. There, positrons are coupled to ma- 
gnetic field lines, which are frozen to the plasma. We distingui- 
shed between advection by star-driven motions and by Galactic 
winds. In the first case, positrons experience the turbulent mo- 
tions of the hot plasma in the winds of massive stars and in super- 
nova remnants. We estimated a diffusion coefficient ~ 3 x 10~^"~^ 
cm- s"', depending on the turbulent velocity and the scale of 
turbulent eddies. In the second case, positrons are transported 
towards high Galactic latitudes by the Galactic wind, over dis- 
tances ~ Vv, Tiife . The Galactic wind velocity ranges from < 1 5 km 
s"' in the Galactic disk to ~200 km s"' a few kpc above the 
Galactic plane. 

If the anisotropy of the magnetic field does not affect the 
scattering off MHD waves in the WIM and HIM, the distances 
traveled by positrons scattered off MHD waves have to be com- 
pared with those resulting from collisions with interstellar mat- 
ter ; the propagation regime leading to the shorter distances is 
the dominant one. Transport by winds will further increase the 
extent of the spatial distributions of positrons. Comparing the 
distances traveled by positrons, we find that scattering off MHD 
waves is more efficient to confine Galactic positrons. For ins- 
tance, assuming quasi-linear diffusion, we estimate that 90% 
of the 1 MeV positrons emitted in the WIM by a point source 
would fill a sphere with radius ~ 20-77 pc before being exclu- 
ded from the resonance condition at an energy ~ 2.5-9.5 keV 
(see Table [1]), depending on the true density and the maximum 
turbulent scale of the WIM. Low-energy positrons would then 
propagate along field lines in the collisional regime before an- 
nihilating. Ultimately, the spatial distribution of these positrons 
would have an extent ranging from ~ 40 to 160 pc, i.e., lower 
than or similar to the typical size of the WIM. This implies that 
all or a large fraction of these positrons would annihilate in the 
WIM, where they are produced. The result is quite different in 
the case of anisotropic turbulence. If we neglect the effect of 
scattering off MHD waves in the WIM, then the spatial distribu- 
tion of positrons would extend between ~ 33 kpc and 100 kpc 
(see Figure [8]), which implies that 1 MeV positrons systemati- 
cally escape from the WIM. 

The effect of the confinement of positrons close to their 
sources has to be taken into account in future investigations 
of their origin. The media in which positrons are likely to be 
produced are generally turbulent (supernova envelopes, super- 
nova remnants, jets, winds). This effect would not only delay the 
time when positrons escape in the ISM, but would also lower 
their "initial" kinetic energy when they are released in the ISM. 
Therefore, such positrons would not propagate too far from their 
birth sites and the spatial distribution of sources would be clo- 
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sely related to the spatial distribution of the annihilation emis- 
sion. Such a scenario is in favor of the hypothesis that LMXBs 
are a dominant source of positrons in the Galactic disk. Indeed, 
since these sources are ex pected to release low-en ergy positrons 
through winds and/or jets (iGuessoum et al 1 120061) . they may ex- 
plain that the longitudinal asymmetry of the annihilation emis- 
sion observed in the Galactic disk fits in with the asymmetry 
observed in the distribution o f LMXBs emitting at high energy 
dWeidenspointner et al. I I2OO8I) . This requires that positrons pro- 
pagate to distances shorter than a few kpc from their sources in 
the Galactic disk, in order to preserve the asymmetry, i.e., that 
their initial kinetic energies are lower than a few hundred keV if 
they propagate in the warm media in the collisional regime. 

To conclude, we showed that the propagation of positrons 
in our Galaxy depends on their initial kinetic energy and on 
the physical conditions of the ISM in which they propagate. 
Therefore, in order to solve the mystery of the origin of Galactic 
positrons, a careful study and modelling of their sources (spa- 
tial distribution in our Galaxy, initial kinetic energy,. . . ) and 
their propagation (ambient medium, mean magnetic field, turbu- 
lence,. . . ) have to be undertaken. Here, we recall that the physi- 
cal parameters of the different ISM phases used in this study are 
characteristic of the Galactic disk. The physical parameters are 
different in the Galactic bulge and particularly close to Sgr A*. 
A detailed calculation for each kind of source will be presented 
in a future paper The results of such calculations will be com- 
pared to the spatial and spectral distributions of the annihilation 
emission measured by SPI. 
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Appendix A: MHD wave cascade damping 

A. 1 . Collisional damping 

Alfven waves 

a. Ionized phases 

With the parameter values listed in Table [1] the proton colli- 
sional mean free-path (Eq.fTSIl is Ap ^ (0.8 - 1.6) x 10'^ cm 
in the hot phase and Ap ^ (1.6 - 5.8) x 10'~ cm in the warm 
ionized phase. Clearly, Linj » /Ip » ^|[J,ax both media. 

In the collisional range Linj ^ ^ ' > -^p, the main dam- 
ping mechanism of Alfven waves is viscous damping (e.g., 
ILazar et al. 20Q3|), which proceeds at a rate 

1 m 



Tvisc ^ - —k^ (sin^ © + 4 cos' ©) , ( A. 1 ) 
2 Pi ^ ^ 



where 



T]\ = 0.3 — r^; + 0.51 



^ce 



(A.2) 



is the (ion + electron) viscosity coefficient dBraginskiil 19651 : 
note that we divided his damping rate, which applies to the 
wave energy, by a factor of 2 so as to obtain the amplitude 
damping rate). When ^ Tg, as will be assumed here, the 
electron contribution to the viscosity coefficient is negligible. 
It then follows from Eqs. ( fT4l i and ( I A. II ) that 



Tvisc Ta < 0.2 ^ 

(airi)2 Va 



1 Vi 1 



In a fully ionized pure-hydrogen gas, the ion r.m.s. velocity, 
Wi = V3 ksTi/nii, is approximately equal to the adiabatic 
sound speed, Cs - yjyPlp. The latter is everywhere smaller 
than or on the order of the Alfven speed, Va - more speci- 
fically, Cs < Va everywhere, except in the hot low-B phase, 
where Cs - (2.7-3.8) Va. Since, in addition, Q.^ Ti » 1, one 
has Fvisc t"a ^ 1 throughout the collisional range (except in 
the limit cos — » 0, which is of little interest here), so that 
the Alfven wave cascade proceeds virtually undamped down 
to Ap. 

b. Atomic phases 

In the weakly ionized (warm and cold) atomic phases, 
Alfven waves are da mped by ion-neutral collisions 
i Kulsrud & Pearcdll969b . Near the injection scale, Alfven 
waves have oj - VA.tot^n ^ Vin,Vni, as can be verified 
with the help of the relations Vin ^ (1.6 x 10"^ cm-' s"')«n 
and Vni - (1.6 X lO"** cm-' s"')«i (see section |2Tt . toge- 
ther with VA,tot = Bl ^ An nip nil, k\\ < 1/Linj, and thus 
VA.tot^ll < (7.1 X 10-i« s-') (b^g/ V«H,cm-4 Again, the re- 
levant parameter values can be found in Table[T] In this low- 
frequency regime, the damping rate due to ion-neutral colli- 
sions is given by 



Tin = 



k^ 

A,tot ni 

2 Vni 



(A.4) 



jFerriere et alj|l988[) . 

Introducing Eq. (IA.4I ) and Eq. ( fT4l i with Va = Va. tot into 
Eq. ( [TtI i then yields for the cutoff parallel wavenumber 



^llcut — 



2Vn 

V, 



= (1.5 X 10""^ cm-') 



V"h,c 



A.tot 



COS 



(A.3) 



(A.5) 

Clearly, this expression is such that VA.tot ~ i^ni ^ Vin, 
which is not quite inside, but not too far from the validity 
Hmit of Eq. (IA.4b . 

A comparison with the right-hand side of Eq. ( fT2l i immedia- 
tely shows that Arycut ^ ^||max, which means that the Alfven 
wave cascade is cut off by ion-neutral collisions way be- 
fore reaching the maximum parallel wavenumber. Moreover, 
since positrons can interact resonantly with Alfven waves 
only over a restricted energy range just above £'k,min, corres- 
ponding to a restricted wavenumber range just below fcytnax 
(see section lTTT i. we may conclude that the Alfven wave cas- 
cade will produce no waves capable of resonantly interacting 
with positrons in the (warm and cold) atomic phases of the 
ISM. 

c. Molecular medium 

In molecular clouds, Alfven waves are again damped by 
ion-neutral collisions. I ncluding the effect of collisi ons of 
gas particles on grains, lElmegreen & Fiebi^ d 1993b accu- 
rately calculated the minimum scale of Alfven waves in a 
molecular cloud of radius R. They found that the maximum 
value of the product Rk is on the order of a few for a 
typical cloud density ~ 10"^ cm"^. Since R typically ranges 
between 0.01 and 10 pc, this leads to a cutoff wavenumber 
^cut ~ (10"'^ - 10"'^) cm"'. Here, too, the smallest scale 
of the Alfven wave cascade is considerably larger than the 
scales at which resonant interactions with positrons occur. 

Fast magnetosonic waves 

a. Ionized phases 

Like Alfven waves, magnetosonic waves with Linj > A:"' > 
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Ap are primarily damped by viscous frictiorQ. However, be- 
cause of their compressible nature, they decay away much 
faster than Alfven waves - by a factor ~ (Qci Ti)-. The cor- 
rect expression of their viscous damping rate reads 



^^k^ sin2 + (rvisc)A 
6 Pi 



(A.6) 



where 



770 = 0.96 «i kBTi Ti + 0.73 «e ksT^ (A.7) 



is the (ion + electron) viscosity coefficient and (rvisc)A de- 
notes the v iscous darnping rate of Alfven waves (given by 
Eq. dAj)) (lBraginskiilll965h . The latter is completely negli- 
gible, except in the limit sin — » 0. 

With both (rvisc)A and the electron contribution to t/o neglec- 
ted, Eqs. (flSl l and ( IA.6I ) combine to give 

Tvisc tf = 0.0533 l}^^ k''^ sin^ © . (A.8) 

As mentioned below Eq. ( |A.3t , v; ^ Cs, so that Vi/Vm^ < 1. 
It then follows that, at the injection scale {k = ^m])' 
Evisc Tf < 1 ■ 

At the transition between the collisional and collisionless 
ranges (k = 



1/2 



[Evisc tf],-. = 0.0533 77^ sin^ , 

Vms \ -^p / 

which can be either smaller or larger than unity, depending 
on the considered ISM phase and on the propagation angle, 
0. 

In the hot phase, [Eyisc tfI^i-i < 1 at all propagation angles. 
This means that the entire fast magnetosonic cascade ma- 
nages to reach the collisionless range beyond A^^, with only 
little or moderate damping. 

In contrast, in the warm ionized phase, [Eyisc tf],-' < ^ oi^ly 
for < 0c, with 0c - 2.60°-4.15°. If all the waves preserve 
their propagation angles across the cascade, then waves with 
< 0c reach the collisionless range, whereas those with 
> 0c are viscous-damped before reaching /tp', and their 
cutoff wavenumber, given by Eq. (fTST l, is a decreasing func- 
tion of : 



^cut = 7.06 



2/3 



-2/3 r-1/3 



L7'" (sin0) 

inj ^ ^ 



-4/3 



(A.9) 



or, numerically, k^^t - ((0.5 - 1.0) x 10""* cm"') (sin ©y"^'^. 
However, the actual situation may not be as clear-cut, as in 
reality, the energy transfer down the cascade is accompanied 
by a randomization of 0, due to both nonlinear interac- 
tions between modes with non-pa rallel wave vectors and 
wandering of magnetic field lines. lYan & LazarianI ( |2004|) 
estimated that the variation in is 50 ~ (kLin^y^^'^. At 
k ~ 10"'"* cm"', this gives 6® ~ 0.02, which is small, but not 
negligible compared to the small values of 0c. Ultimately, 
we may not rule out the possibility that even waves with 
< 0c are cut off by viscous damping before reaching A^^ 



We checked that Joule damping leads to a cutoff wavenumber se- 
veral orders of magnitude greater than zip'. It only dominates viscous 
damping at propagation angles — > 0. 



b. Atomic phases 

Here, fast magnetosonic waves are damped by ion-neutral 
colUsions, at a rate 



2Vni 



/(©), 



(A. 10) 



with 



/(©) = 



1 



1 + 



n,ot + c^-2yi.„.c2 cos 20 



A.tot 



(lEerriereetalJll988h . 

Substitution of Eqs. dA.lOl l and ( fTSl l into Eq. (fTSb directly 
gives for the cutoff wavenumber 



kc\\\ — 



2 Vni 

V A.tot 



-1/2 



/(0)j 



2/3 



(A. 11) 



Since VA.tot > Cs in both atomic phases, we may, to the order 
of the present approximation, let V^s - VA.tot and /(0) - 1 
in the above equation, whereupon we obtain 



^cut — 



2 Vn 



' A.tot 



2/3 



-1/3 



inj 



(0.9 X 10"'Vm"') 



2/3 



(A. 12) 



The fast magnetosonic wave cascade is cut off by ion-neutral 
collisions way before reaching the small scales at which po- 
sitrons can be in resonant interaction. 

c. Molecular medium 

Our conclusion is similar to that reached for Alfven waves. 



A.2. Collisionless damping 



In section lATI we saw that, in the ionized phases of the ISM, 
the Alfven wave cascade experiences negligible collisional dam- 
ping, which enables it to make it all the way down to the colU- 
sionless range ^"' < Ap. The fast magnetosonic wave cascade, in 
contrast, experiences significant collisional (viscous) damping. 
In the hot ionized phase, the collisional range Linj > A:"' > /Ip is 
sufficiently narrow that the fast magnetosonic cascade neverthe- 
less reaches the collisionless range only partially attenuated. But 
in the warm ionized phase, where the collisional range spans 
roughly eight decades, the fast magnetosonic cascade comple- 
tely (or almost completely) decays away before reaching the col- 
lisionless range. 

We now examine collisionless damping in the cases of inter- 
est, namely, for the Alfven cascade in the hot and warm ionized 
phases and for the fast magnetosonic cascade in the hot ionized 
phaselfl In all cases, the dominant collisionless damping mecha- 
nism is linear Landau damping. 

Alfven waves 



^ We will, however, keep in mind the possibility that, in the warm io- 
nized phase, waves with small propagation angles might enter the colli- 
sionless regime. 
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For Alfven waves, a n approximate e xpression of the linear 
Landau damping rate is ( iGinzburdl 19611) : 



LD 



71 , COS sin^ 

sin- + 3 — cos"* 



+ (sin^ + 4 cos^ 0) exp | 



y2 ^ 



(A.13) 



where = yJk^TJmi and v'^ - yfk^Tjml are the ion and elec- 
tron thermal speeds (which differ from the ion and electron r.m.s. 
velocities, v\ and u^, by a factor V3) and k/^ - Q^i/^A is the 
inverse ion inertial length. For a pure-hydrogen plasma, kp^ is 
nothing else than the maximum parallel wavenumber, ^ymax (see 
Eq. (0)). 

Multiplying Eq. ( IA.13I ) by Eq. ( fT4l l immediately leads to 



Fld ta 



n v' k- 



sin2 



l<l . . k^ . 

sm' + 3 -r- cos 

^1 



j2 + (sin^ + 4 cos^ 0) exp 



Vi 



2v' 



i / 



(A. 14) 



The cutoff of the Alfven cascade occurs when Fld ta = 1 (see 
Eq. (fTTIi). provided that this relation admit a real solution for k. 

In the hot low-B and warm ionized phases, where Va ^ f[, 
the first term inside the square brackets on the right-hand side of 
Eq. ( IA.14b is negligible, and the relation Fld ta = 1 admits one 
solution, given by 

, ISV'' /VaV'' 
^cut ^ - — exp 



(sin^ + 4 cos^ ©) ^'^ k 



(A.15) 

(except in the limit sin — > 0). Ignoring the weak 0-dependence 
of kcut (which entails only a variation by a factor of 2), we find 
kcnt - (0.14 - 0.18) in the hot low-B phase and k^ut - (0.65 - 
3.20) ^A in the warm ionized phase. Hence, in these two media, 
the Alfven cascade is cut off by Landau damping at a scale close 
to the proton inertial length. 

In the hot high-Z? phase, where Va exceeds i/| by a fac- 
tor ^ 4.8 - 6.8, the exponential factor on the right-hand side 
of Eq. (IA.I4I 1 becomes negligibly small, with the result that 
Eld Ta < 1 at all k (except in the limit cos © — > 0, but in this 
limit Alfven waves are viscous-damped before reaching the col- 
lisionless range ; see Eq. (IA.3l l). If taken at face value, this result 
would lead to the erroneous conclusion that the Alfven cascade 
is not cut off by Landau damping. 

The truth is that our expressions for the frequency, transfer 
time and damping rate of Alfven waves are all strictly valid only 
in the MHD regime, so that Eq. ( IA.14b actually breaks down 
at wavenumbers approaching the inverse proton inertial length, 
^A- The only conclusion that can be drawn from Eq. ( |A.14| i is 
that the cutoff of the Alfven cascade by Landau damping occurs 
roughly at the inverse proton inertial length, ^a- 

Fast magnetosonic waves 

For fast m agnetosonic wav es, linear Landau damping pro- 
ceeds at a rate (IGinzburdl 96 ll) 



Frn = 



7T , , Sin 

v\ k 

cos © 



-J +5 exp 



Vi 



2yp cos^©^ 



which, combined with Eq. ( fTsT l. yields 



1/2 



sin^0 
cos 



-J- + 5 exp 

v' 



Vi 



2 v'~ cos^ © 
' (A. 17) 

The cutoff wavenumber, at which Fld tf = 1 (see Eq. (fTSb). is 
given by 



8 



COS- 

sin"*© 



— + 5 exp 



y2 



2 v'"- cos2 © j 

1 ' 



inj 



> i/j, the exponen- 
( IA.18b . and the result 

lO^Lr^] cos2©/sin'*© ^ 



(A. 16) 



(A.18) 

Here, the cutoff depends strongly on ©, and it diverges at paral- 
lel propagation (just like the viscous cutoff in the warm ionized 
phase ; see Eq. ( IA.9I )). Once again, due to the possible rando- 
mization of (see section IATTT i. this mathematical divergence 
should not be a concern, unless the randomization process is in- 
efficient. 

The only interstellar phase that needs to be considered here 
is the hot ionized phase. In the low-B case, Va < u\ and the 
expression within square brackets in Eq. ( IA.18b is dominated by 
the exponential term (except in the limit cos © 0). Eq. ( IA.18I) 
then reduces to fccut ^ (0.46 - 0.61) Lr^j cos- 0/ sin"* 0. Despite 
the strong ©-dependence of k^y^^, we may conclude that the fast 
magnetosonic cascade is severely affected by Landau damping 
as soon as it enters the collisionless range and that it globally 
decays away shortly below A^. 

In the high-B case, Va 
tial term drops out from Eq 
reads A;cut - (1-2 - 2.3) x 

((4.0 - 7.5) X 10-'^ cm-i) cos- 0/ sin"* ©. Here, Landau 
damping is less severe, which allows the fast magnetosonic 
cascade to globally proceed over roughly three decades below 
/Ip, before vanishing. 
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